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Abstract 

I present a calculation of QCD radiative corrections to charged Higgs production via 



o 
o 



the process bg — > tH~. I show that the cross section is dominated by soft-gluon cor- 
, rections, which are computed through next-to-next-to-leading order. Results for charged 

q ' Higgs production at the LHC are presented, including the dependence of the cross section 

on the charged Higgs mass, the top quark mass, the factorization and renormalization 
scales, and tan (5. 
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1 Introduction 

> 

One of the main goals of the current particle physics program is the discovery of the Higgs boson. 
The Minimal Supersymmetric Standard Model (MSSM) introduces charged Higgs bosons in 
addition to the neutral Higgs. A future discovery of a charged Higgs would thus be a sure sign 
of new physics beyond the Standard Model. 

In the MSSM there are two Higgs doublets, one giving mass to the up-type fermions and 
the other to the down-type fermions. The ratio of the vacuum expectation values, v 2, V\ for the 
two doublets is tan/5 = V2/V1. Among the extra Higgs particles in the MSSM are two charged 
Qh! Higgs bosons, H + and H . 

The Large Hadron Collider (LHC) at CERN has a good potential for discovery of a charged 
Higgs boson. A promising channel is associated production with a top quark via bottom-gluon 
fusion, bg — > tH~ [1-17]. In this paper we focus on H~ production, but we note that the 
cross sections for H + production, via the related process bg — > iH + , are identical. The 
complete next-to-leading order (NLO) QCD corrections to the process bg — > tH~ have been 
recently derived in Refs. [T2] ITS] IT3j . These corrections were shown to stabilize the cross 
section with respect to changes in factorization and renormalization scales. The SUSY-QCD 
NLO corrections were also calculated in [T51 I14j. 

The NLO QCD corrections were shown to be substantial, up to 85% enhancement of the 
lowest order cross section [T2j. The SUSY-QCD corrections are comparatively small, though 
non-negligible, and their precise value depends on several parameters of the MSSM ^H]- Since 
the NLO QCD corrections are large it is important to consider whether even higher-order 
corrections may make a significant contribution. In this paper I show that the NLO corrections 
are dominated by near-threshold soft-gluon emission and I calculate the contribution from 
next-to-next-to-leading order (NNLO) soft-gluon corrections, which are sizable. 

The charged Higgs is expected to be quite massive, so its production at current colliders will 
be a near-threshold process. In such processes soft gluon emission is expected to dominate the 
radiative corrections. This has in fact been shown by now for a large number of processes in- 
cluding top, bottom, and charm quark production [TH], H^-boson production [19], direct photon 
production [20], jet production [21], and flavor-changing-neutral-current single-top production 
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|22| . Near threshold for the production of a specified final state there is limited energy available 
for the production of any additional radiation; hence the emitted gluons are soft and they man- 
ifest themselves in logarithmically enhanced terms that numerically dominate the cross section. 
The structure of these threshold contributions follows from general considerations of the factor- 
ization properties of hard-scattering cross sections [23 1211 1213 ■ Renormalization properties of 
the factorized pieces of a cross section result in formal resummations, which provide the form 
of the soft-gluon corrections to all orders in the strong coupling a s . For further details and 
reviews see Refs. [221 123 EES • 



2 NNLO soft-gluon corrections 

In this section we derive the analytical form of the soft-gluon corrections through next-to-next- 
to-leading order for charged Higgs production with a top quark in hadronic collisions. For the 
process b(pb) + g(j>g) ~~ > t(Pt) + H~(Ph~), we define the kinematical invariants s = (pb + p g ) 2 , 
t = {pb — Pt) 2 , u = (p g — p t ) 2 , and S4 = s + t + u — m 2 — where m H - is the charged 

Higgs mass and m t is the top quark mass. Note that we ignore the mass of the 6-quark in 
the kinematics. Near threshold, i.e. when we have just enough partonic energy to produce 
the tH~ final state, S4 — > 0. The threshold corrections then take the form of logarithmic 
plus distributions, [(ln'^/m 2 ^)/^)]^ where I < 2n — 1 for the n-th order QCD corrections. 
These plus distributions are defined through their integral with any smooth function, such 
as parton distributions, giving a finite result. The leading logarithms (LL) are those with 
I = 2n — 1 while the next-to-leading logarithms (NLL) are those with I = 2n — 2. In this 
paper we calculate NLO and NNLO soft-gluon threshold corrections at NLL accuracy, i.e. at 
each order including both leading and next-to-leading logarithms. We denote them as NLO- 
NLL and NNLO-NLL, respectively. Thus, at NLO we include [ln(s4 : /m 2 H J)/s4\ + (LL) and 
[l/s 4 ] + (NLL) terms. Although we do not calculate the full virtual S(s4) terms, we include 
those 5(54) terms that involve the factorization and renormalization scales, denoted by [if and 
fin respectively. At NNLO, we include [ln 3 (s4/m^-_)/s4] + (LL) and [In 2 (s4/m 2 H -) / 34]+ (NLL) 
terms. We also include some [In(s4/m 2 ? _)/s4] + and [l/s 4 ] + terms that involve the factorization 
and renormalization scales; and some constants which arise from the inversion from moment 
space, where the resummation is performed, back to momentum space. For details of this 
approach see Refs. [T%1 HF]. 

The differential Born cross section is d 2 a b B g ^ tH - j \dt du) = F h B g _ ttH -8{s^) where 

B rraa s (ml tan 2 (3 + m 2 cot 2 (3) j s + t — m 2 H - 



b9 ^ tH ~ ' I2s 2 m 2 w sm 2 9 w \ 2s 

m 2 (u — m 2 H J) + m 2 H - (t — mf) + s(u — mf) m 2 (u — m 2 H - — s/2) + su/2 



s{u — m 2 ) {u — wq) 



2\2 



,(2-1) 



where a = e 2 /(4n), a s is the strong coupling, and we have kept the 6-quark mass, mj,, non-zero 
only in the m 2 tan 2 (3 term. We use consistently the running masses for the top and bottom 
quarks j2Ej, corresponding to pole masses of 175 GeV and 4.8 GeV, respectively. 

We next proceed with the calculation of the NLO and NNLO soft-gluon corrections at NLL 
accuracy. In our derivation of these corrections we follow the general techniques and master 
formulas presented in Ref. 
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The NLO soft-gluon corrections for the process bg — > tH are 



d 2 a {1) 
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- 2(CV + CU), where C F = (iV 2 - 1)/(2N C ) and C A = iV c with N c = 3 the 
number of colors, and 
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where fi F is the factorization scale, and we have defined j 1 " 9 ^^ as the scale-independent part 
of c b 2 9 ^ tH . The term Rel 1 '^ denotes the real part of the one-loop soft anomalous dimension, 
which describes noncollinear soft-gluon emission , modulo some gauge-dependent terms that 
cancel out in the cross section. A one-loop calculation gives 
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where (Ir is the renormalization scale and A) = (11C A — 2n/)/3 is the lowest-order j3 function, 



with rif the number of light quark flavors. Note that c 



bg^tH- 

l 



represents the scale-dependent part 



of the S(s^) corrections. We do not calculate the full virtual corrections here. Our calculation of 
the NLO soft-gluon corrections includes the full leading and next-to-leading logarithms (NLL) 
and is thus a NLO-NLL calculation. 

We next calculate the NNLO soft-gluon corrections for bg — > tH~: 
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where (2 = 7r 2 /6 and (3 = 1.2020569.... We note that only the leading and next-to-leading 
logarithms are complete. Hence this is a NNLO-NLL calculation. Consistent with a NLL 
calculation we have also kept all logarithms of the factorization and renormalization scales in 
the [In(s4/m^_)/s4] + terms, and squares of logarithms involving the scales in the [1/54]+ terms, 
as well as (2 and (3 terms that arise in the calculation of the soft corrections when inverting 
from moments back to momentum space [T%1 125j. 

In principle one can obtain the form of the soft radiative corrections at any order in a s and 
indeed resum them to all orders. However in practice such resummed cross sections depend on 
a prescription to avoid the infrared singularity and ambiguities from prescription dependence 
can actually be larger than contributions from terms beyond NNLO ^H]- Hence we here give 
results to NNLO as has been done for many other processes [EH CHI 033 1201 I2H 1221 122] ■ 

We now convolute the partonic cross sections with parton distribution functions to obtain 
the hadronic cross section in pp collisions at the LHC For the hadronic cross section p(pa) + 
p(Pb) -» t(p t ) + H~(p H -) we define S = (p A + Pb) 2 , T = (p A - p t ) 2 , and U = (p B - Pt) 2 , 
and note that pb = x A p A , p g = XbPb, where x denotes the momentum fraction of the hadron 
carried by the parton. The hadronic cross section is then given by 
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where 



x A 



s 4 — m 2 + m 2 H ^ — Xb(U — m 2 
XbS + T — mf 



X8) 



Trna* = - (1 / 2) (S - m* - w^-) ± (1/2) J (S + mf - mjj.) 2 - AmjS, and <j)(x) are the parton 
distributions. 



3 Charged Higgs production at the LHC 

We now turn our attention to detailed numerical results for charged Higgs production at the 
LHC. In Figure H we plot the cross section versus charged Higgs mass for pp collisions at the 
LHC with \/~S = 14 TeV. Here and throughout this paper we use the MRST2002 approximate 
NNLO parton distributions |2Z| with the respective three-loop evaluation of a s . We set the 
factorization scale equal to the renormalization scale and denote this common scale by \x. We 
show results for the Born, NLO-NLL, and NNLO-NLL cross sections, all with a choice of scale 
fi = mil-. We use the same NNLO parton densities and couplings in all the results, so that we 
can concentrate on the effects of the soft-gluon corrections. We choose a value tan/3 = 30. It is 
straightforward to get the results for any other value of tan /3, since the only dependence on /3 
in our equations is in the factor ml tan 2 j3 + m\ cot 2 (3 appearing in the Born term, Eq. 12.11 The 
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Figure 1: The total cross section for charged Higgs production at the LHC. 



cross sections span over two orders of magnitude in the mass range shown, 200 GeV < van- < 
1000 GeV. The NLO and NNLO threshold corrections are positive and provide a significant 
enhancement to the lowest-order result. We note that the cross sections for the related process 
bg — > iH + are exactly the same. 

The relative size of the corrections is better shown in Figure El where we plot the i^-factors, 
i.e. ratios of cross sections at various orders. The NLO-NLL / Born curve shows that the NLO 
threshold corrections enhance the Born cross section by approximately 25% to 50% depending 
on the mass of the charged Higgs. As expected the corrections increase for higher charged 
Higgs masses since then we get closer to threshold. The NNLO-NLL / Born curve shows that 
if we include the NNLO threshold corrections we get an enhancement over the Born result of 
approximately 35% to 70% in the range of masses shown. Again the enhancement increases 
with charged Higgs mass, as expected. Finally, the NNLO-NLL / NLO-NLL curve shows clearly 
the further enhancement over NLO that the NNLO threshold corrections provide. This curve 
is simply the ratio of the other two curves and varies between 7% and 14%. 

We now want to compare our NLO-NLL results with the exact results that have been derived 
in |12[ [T3] . We note that different choices of factorization/renormalization scales were used in 
those references. In Ref. ^2] the reference scale chosen was van- + m t while in Ref. [T3| it was 
{mu- + n2 t )/2. In this paper we choose van-- This is the natural choice in our approach since 
we are considering logarithms of s^jva 2 H -. Of course any choice of scale is theoretically possible 
and a cross section known to all orders does not depend on the scale. However a finite-order 
cross section does depend on the scale, though the dependence decreases as we move from 
Born to NLO to NNLO and so on. The work in ^21 E21 indeed showed a reduction of scale 
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Figure 2: The i^-factors for charged Higgs production at the LHC. 



dependence when the NLO corrections are added relative to the Born cross section. In fact, as 
we will see below, the NNLO threshold corrections further decrease the scale dependence, thus 
resulting in more stable predictions. 

Before comparing our results to the exact NLO cross section, we can check the effect of 
choosing the scales used in Refs. (TJHinj. In Figure El we plot the ratios of the cross sections 
with choice of reference scale fj, = ran- + rn t (bold lines) and \i = (m H - + m t )/2 (thin lines) 
over the cross section with scale /x = m H - . We see that indeed there is a considerable variation 
at lowest order, but this progressively diminishes at NLO and NNLO. In fact at NNLO there 
is hardly any difference between the two lines at large values of the charged Higgs mass. Thus 
we see the stabilization of the cross section versus scale variation when higher-order corrections 
are included. We will say more regarding this important point and show more plots below. 

We now compare the NLO soft-gluon results with the exact NLO cross section. In Figure 
El we compare the NLO-NLL cross section with the exact NLO cross section of reference |12j . 
To make the comparison, the NLO-NLL result is calculated here for \x = tuh- +mt since that's 
the scale chosen in [TJ] and also using a two-loop a s . Also, to remove discrepancies arising 
from different choices of parton distribution functions, we plot iT-factors. The NLO-exact / 
LO curve is taken from Ref. ^2] by dividing curve 1 by curve 2 in Figure 6 of that reference 
(to account for the different definition of i^-factor used there). The fact that the NLO-NLL / 
NLO-exact curve is very close to 1 (only a few percent difference) shows that the NLO-NLL 
cross section is a remarkably good approximation to the exact NLO result. As noted before, 
we might have expected this on theoretical grounds since this is near-threshold production, 
and also from prior experience with many other near-threshold hard-scattering cross sections 
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bg->tH at LHC S =14TeV n=m H +m t , (m H -+m t )/2 
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Figure 3: The ratio of cross sections at various scales for charged Higgs production at the LHC. 
The bold lines are with /i = + m t ; the rest are with /i = (ra#- + m t )/2. 
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Figure 4: The ratio of exact and approximate NLO cross sections for charged Higgs production 
at the LHC. 
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bg --> tH" at LHC S =14TeV tan|3=30 m -=500 GeV 
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Figure 5: The scale dependence for production of a charged Higgs with mass van- — 500 GeV 
at the LHC. 
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In Figure El we plot the scale dependence of the cross section for a fixed charged Higgs mass 
m H - = 500 GeV and tan/? = 30. We plot a large range in scale, 0.1 < fi/m H - < 10, and see 
indeed that the threshold corrections greatly decrease the scale dependence of the cross section. 
The NNLO-NLL curve is relatively flat. 

In Figure El we plot the the cross section as a function of charged Higgs mass with two 
different choices of scale, \i = van- /2 and 2m#-. We see that the variation with scale of the 
Born cross section is quite large. The variation at NLO-NLL is smaller, and at NNLO-NLL it is 
very small. In fact the two NNLO-NLL curves are on top of each other for most of the range in 
run- ■ Hence, the scale dependence of the cross section is drastically reduced when higher-order 
corrections are included. This is as expected from and is consistent with the reduced scale 
dependence shown in Figures 01 and El 

In Figure [7| we plot the dependence of the cross section on tan (3, over the range 1 < 
tan/3 < 50, for fixed charged Higgs mass and scale \i = m H - = 500 GeV. The cross section 
is at a minimum near tan/5 = 8. We note that the tan/? dependence arises in the factor 
va\ tan 2 (3 + m\ cot 2 (3 in the Born term and thus the shape of the curves at higher orders is 
similar. The dependence on tan/3 is quite large, spanning nearly two orders of magnitude in 
the range shown. 

In Figure |H1 we plot the dependence of the cross section on the top quark mass for fixed 
charged Higgs mass and scale /i = m H - = 500 GeV and tan (3 = 30. We see that the dependence 
is not very strong so that the present experimental uncertainties on the top quark mass do not 
play a dominant role in the total uncertainty of the charged Higgs production cross section. As 
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Figure 6: The scale dependence for charged Higgs production at the LHC The bold lines are 
with jj, = 2rriH-', the rest are with \i = mu-jl. 



bg --> tH" at LHC S =14 TeV (i=m H =500 GeV 

1 r — i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



Born 

NLO-NLL 

-I NNLO-NLL 




tan (3 

Figure 7: The tan (3 dependence for charged Higgs production at the LHC. 
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Figure 8: The top quark mass dependence for charged Higgs production at the LHC. 
the top quark mass gets more precisely known, this dependence will diminish further. 

4 Conclusion 

The process bg — > tH~ offers a promising possibility for discovering a charged Higgs boson. 
Charged Higgs production at the LHC receives important contributions from the threshold 
region. The NLO corrections to the process bg — > tH~ are quite large. We have seen that the full 
NLO cross section is very well approximated by the NLO-NLL soft-gluon result, to within a few 
percent. The NNLO soft-gluon threshold corrections to charged Higgs production are important 
and further stabilize the cross section versus changes in factorization and renormalization scales. 
The dependence on tan (3 and on the charged Higgs mass are quite large while the dependence 
on the top quark mass is milder. 
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